Softening induced instability of a stretched cohesive granular layer 
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We report on a cellular pattern which spontaneously forms at the surface of a thin layer of a 
cohesive granular material submitted to in-plane stretching. We present a simple model in which 
the mechanism responsible of the instability is the "strain softening" exhibited by humid granular 
materials above a typical strain. Our analysis indicates that such type of instability should be 
observed in any system presenting a negative stress sensitivity to strain perturbations. 
PACS: 89.75.Kd; 83.60.Uv; 45.70.Qj. 



Introduction - Adding even minute amounts of liq- 
uid can change dramatically the mechanical properties of 
sand. During the building of sand castles, one observes 
a transformation from a fluid-like to a sticky and de- 
formable material with increasing water content. Indeed, 
at very low water content, the formation of partially- 
developed capillary-bridges leads to a fast increase of 
tensile strength whereas, for large enough fluid content, 
tensile strength is nearly constant Cluster formation 
was identified as the main mechanism responsible of such 
a behavior 

One unexplored but important feature of the mechan- 
ical response of cohesive sand, arising in a wide range of 
fluid content, is that the tensile strength decreases when 
the imposed strain is increased. Such a "strain softening" 
is due both to a decrease of the associated adhesion force 
when a single bridge is elongated [3] and to a decrease 
in the overall number of bridges which collapse when ex- 
cessively stretched [4J. We mention that the softening 
behavior is only observed above a critical strain which 
vanishes for sufficient stiff grains Q. Indeed, the critical 
strain is associated with the initial compression of the 
grains induced by the suction force due to the capillary 
bridges. 

In the present Letter, we report that, in a stretched 
layer of cohesive grains, "strain softening" induces a me- 
chanical instability in which the strain field is modulated 
in space. Our analysis indicates that such type of in- 
stability should be observed in any system presenting a 
negative stress sensitivity to strain perturbations. 

Experimental setup and protocol - The experiment 
consists in imposing an in-plane deformation at the base 
of a thin layer of a cohesive granular material. To do so, 
the grains are initially spread onto an elastic membrane 
to which the deformation is imposed (Fig. [T]). The mem- 
brane, a thin Latex band (thickness 0.5 mm, width 5 cm 
and total length 20 cm) is maintained at its two ends by 
a moveable U-shaped frame and leans, in its central part, 
on a steady rectangular table (length 10 cm, width 6 cm). 
By displacing the frame downwards, the band, which re- 



mains in the same horizontal plane above the table, ex- 
tends along its length and narrows in the perpendicular 
direction (Note that, due to the contact with the table, 
the membrane does not wrinkle). The top plate of the 
table is made of polytetrafluoroethylene (PTFE), which 
limits the friction and insures a reproducible deformation 
for a given frame displacement. 
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FIG. 1. Sketch of the experimental setup - Inset: Top view 
of the whole pattern and definition of the angle a. 

The sample is prepared first by pouring grains onto 
the membrane. The surface of the material is then lev- 
eled by means of a cylindrical rod guided by lateral spac- 
ers, which achieves a well-defined thickness h (from 1 
to 5 mm, to within 0.1 mm). The granular material con- 
sists of spherical glass-beads (USF Matrasur, sodosilicate 
glass). We shall report results obtained for various sam- 
ples in a large range of bead diameters d (0-45, 45-90, 
100-125 /xm). In order to tune the cohesion, the experi- 
mental device is placed in a chamber in which the atmo- 
sphere is equilibrated with a saturated salt solution. The 
cohesion is, in addition, accounted for by measuring the 
angle of avalanche 9 a @,0| m the same experimental con- 
ditions (A granular layer, same material and thickness, 
is prepared onto a rough surface which can be tilted). 
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The free surface of the sample is imaged from above by 
means of a digital camera (Konica-Minolta, A200). Two 
linear light-sources (home-made arrays of LEDs, Fig. [T]) 
placed at the two ends of the elastic band, about 20 cm 
away from the sample, a few centimeters above the table 
plane, provide a good contrast when the upper surface of 
the material is deformed. In addition, in order to assess 
quantitatively the vertical displacement of the free sur- 
face, we cast with an angle of 30 deg, a light sheet onto 
the sample. In this configuration, the horizontal displace- 
ment of the bright line is proportional to the local vertical 
displacement of the free surface. 

Results - When the membrane is stretched by mov- 
ing the frame downwards, one observes, provided that 
the grains are small enough and/or the relative humid- 
ity, Rh , large enough, the growth of a complex pattern 
at the free surface of the granular layer (Fig. [2]). 
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FIG. 2. Surface structure for distinct layer thickness and 
particles size at constant stretching - Lower left panel : a 
highly cohesive layer slides on the membrane (Labels: Top, 
h; Right, d; 6 ~ 0.2) 

Domains, made of stripes having a rather well-defined 
width and making an angle ±ev with the stretching direc- 
tion (x-axis, Fig. [I]), nucleate and grow. We note that the 
phenomenon consists of the continuous growth of a free- 
surface undulation, which, depending on the experimen- 
tal conditions, can eventually lead to the fracture of the 
granular layer when the stretching is further increased. 
Note that the undulation first appears at the edges of 
the latex band where it makes a 90 deg angle with the 
stretching direction (Fig.[TJ inset), which indicates, first, 
that the uniaxial stretching can be responsible, alone, for 
the instability, and, second, that the compression along 
the y-axis plays a role in the angle a measured in the 
central part of the sample. We can also deduce from the 
observation of domains, disconnected from the edges, ex- 
hibiting a given orientation (±a) that the overall pattern 



does not result from the growth of the instability occur- 
ring at the lateral boundaries. 

Let us first discuss the angle a. We denote 9 = 
u xx , the strain in the x direction (u denotes the dis- 
placement field) and report a measured for a given 9 
(Fig. [3^). We observe that, to within the experimen- 
tal accuracy (±3 deg on one sample), a neither depends 
on h, d or cohesion (accounted for by 9 a ) and we get 
a = (51.5 ± 0.5) deg for 9 ~ 0.21. In contrast, a de- 
pends on 9. Indeed, the latex being almost incompress- 
ible, u yy ~ —9/2 and the stripes are expected to rotate 
according to tan(a) = 2 fi+e) taua(ao). Reporting tan(a) 
as a function of 9 (Fig. [3b), we get aa ~ 58 deg, the 
value of a in the limit 9 — > 0. One can account for 
the experimental angle ckq by considering the Mohr cri- 
terium Q . On the one hand, the stretching tends to pull 
the grains apart, so that the normal stress along the x- 
axis, a xx < a s , where a s denotes the tensile stress. On 
the other hand, the compression pushes the grains one 
against another and one can guess that the associated 
stress involves the solid contacts between the grains and, 
consequently, that \<J yy \ S> |o"xx|- Thus, one can consider 
that the granular material is subjected to a pure compres- 
sive stress along the y-axis. With this assumption, the 
Mohr criterium leads to «o = f + % where tan ($) = \i, 
the static friction coefficient ||. The experimental value 
of ao corresponds to $ ~ 26 deg associated with /j, ~ 0.5, 
a reasonable value of the static friction coefficient for a 
packing of glass spheres @ . 
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FIG. 3. (a) Angle a vs. thickness h for 9 ~ 0.21 - The angle 
a neither depends on the thickness h of the granular layer nor 
on the grain size d and cohesion [Symbol : d (9 a )]- (b) Angle 
a vs. elongation 9 - The stretching induces a rotation of the 
stripes such that tan(a) = (2 — 0)/[2(l + 9)] tan(ao). From 
the extrapolation of the data (grey line), one assesses the 
angle ao ~ 58 deg (d =0-45 /im, h — 2 mm and 9 a — 47 deg). 

We measured the typical width of the stripes and am- 
plitude of the vertical deformation of the layer. We ob- 
served that similar results are obtained when a pure uni- 
axial deformation of the membrane is imposed. Thus, 
seeking for simplicity, we shall report measurements as- 
sociated with the deformation of the layer in the simplest 
geometry: we limit the experiment to 1-cm-wide mem- 
brane stripe (length 5 cm), limited by two thin metal 
plates glued onto the membrane (Fig. @Ji)- 

When the membrane is stretched, the central region, 
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FIG. 4. (a) Sketch of the experimental configuration for uni- 
axial strain, (b) Photograph of the fracture pattern - Frac- 
tures are perpendicular to the strain axis. The white line is 
the vertical displacement u z (x,h) of the surface, (c) Verti- 
cal displacement of the free surface u z (x,h) - The depth at 
the fracture increases whereas the free surface is observed to 
rise, and thus the granular layer to dilate, between the frac- 
tures when 9 is increased (A darker line indicates a larger 9) . 
(d) Wavelength A (diamonds) and amplitude (squares) of the 
vertical undulation u z (x,h) vs. strain 9 - The amplitude 
increases continuously with 9 whereas, due to the simple ad- 
vection of the pattern, A obeys Ao(l + 9) (d — 100-125 /im, 
h — 2 mm and 9 a = 37 deg). 



away from the edges, is subjected to a pure uniaxial strain 
and one then observes the formation of stripes perpen- 
dicular to the elongation axis (Fig.0p). The deformation 
of the bright line casted onto the surface (Fig. [T]) makes 
it possible to obtain the vertical profile of free surface 
along a line and thus to quantify the growth of the pat- 
tern (Fig.Ht). We observe that, obviously, the fracture is 
marked by a sharp decrease of the surface altitude but, 
interestingly, that the height of the free surface in the 
domains between the fractures increases and, thus, that 
the in-plane stretching leads to a thickening of the layer. 

At this point, it is interesting to consider the depen- 
dency of the typical wavelength, A, of the fracture pat- 
tern on the experimental parameters that are h, d and 
Rh- Note first that A is not strictly selected and that 
we observe a large scatter of the stripes width. In spite 
of the scatter, we observe that A oc h (Fig. [SJ as long 
as h does not exceeds about 3 mm (If h is too large, 
the bands between the fracture are likely to split in two, 
leading to a smaller average A). For a given Rh and h, A 
is almost independent of d (Fig. [5]) whereas one observes 
a significant increase of A when cohesion is increased by 
increasing Rh (Inset: Fig. [SJ- 

Theoretical analysis - The instability requires cohe- 
sion. Due to the nature of the interaction between the 
grains, the adhesion force decreases when the material is 
stretched and, thus, grains are pulled apart 0,0]. 
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FIG. 5. Wavelength A vs. thickness h - Full squares: A 
is proportional to h. For d = 45-90 /im, A ~ 0.6 h; Open 
squares: A does not significantly depend on d (9 ~ 0.3 and 
Rh ~ 35%). For clarity, error bars (20% in A, 0.1 mm in 
h) are not shown. - Inset: A increases with Rh (9 ~ 0.3, 
d — 100-125 /im and h = 2 mm). 



Qualitatively, when the membrane is elongated, an ho- 
mogeneous stretching of the material is imposed in the 
base plane. However, due to the "strain softening" , in 
response to the overall stretch, the system tends spon- 
taneously to modulate the deformation: regions of large 
deformation are associated with a smaller modulus and 
regions of large modulus are associated with a smaller de- 
formation, which results in an overall decrease of the en- 
ergetic cost. In turn, the modulation induces a shear de- 
formation which is associated to an energetic cost. Thus, 
the wavelength is governed by the balance of the gain 
associated with the modulation of the horizontal strain 
and of the loss associated with the induced shear. We 
point out here, that the shear (the relative motion of the 
top and bottom surfaces) is consistent with the dilation 
of the layer in the vertical direction (Fig. 

In order to account, at least qualitatively, for the ex- 
perimental observations let us first note that the normal 
stress along the decreases linearly with the 

uniaxial strain u XXl according to <r xx = cr s (l — u xx /9 m ) 
when the material is stretched (u xx > 0) [ajj. The rela- 
tion is no longer valid for u xx > 9 m , when the elongation 
is large enough for the bridges to collapse and, thus, the 
material to break apart. Thus, 8 m is of the order of the 
typical size of the bridge S divided by the grains diameter 
d whereas <j s denotes the tensile stress previous to defor- 
mation. For the sake of simplicity, the contribution of the 
shear shall be accounted for by a simple shear modulus G 
which value shall be discussed later. In this framework, 
the shear stress o xz — Gu xz and, accordingly, the energy 
per unit volume 



/ u \ 1 

E = <7. [u xx — + -Gu 

\ 26U J 2 



2 

xz ' 



(1) 



We thus assume that the dilation in the vertical direc- 
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tion u zz does not contribute to any additional energetic 
cost: The grains are assumed to remain in contact along 
the vertical and no significant stretching of the capillary 
bridges occurs in this direction. 

We consider now a sinusoidal perturbation of the dis- 
placement such that u x = Ox + f(z) sin (kx) in the hori- 
zontal plane. In order to obtain the associated displace- 
ment in the vertical direction u z , seeking for simplicity, 
we will further assume that the deformation of the ma- 
terial does not induce any dilation so that u xx + u zz = 0. 
Writing the relation imposed by the mechanical equilib- 
rium in the horizontal plane + ^*j%r = [13 and the 
conditions that u x and u z do not depend on x at the sub- 
strate plane (z = 0), we get u x = Qx + auj sin {ujz) sin (kx) 
and u z = —9z — ak[l — cos (cuz)] cos (kx) where a is an am- 
plitude and lo 2 = k 2 (l+C) with £ = 2a s /(G0 m ). Writing 
that the shear stress a xz vanishes at the free surface and, 
thus, that u xz — for z — h, we get 1 + £ cos (u)h) = 0. 
The wavelength A = 2n/k is thus found to be propor- 
tional to h, independent of 9, according to 



A = 2tt 



arccos (—!/£) 



(2) 



provided that £ > 1. One can check that in the accessible 
range of £, the energy E is a decreasing function of the 
amplitude a, whatever the strain 9. Thus, the layer is 
always unstable provided that the decrease in the tensile 
stress is large enough compared to the shear cost, i.e. 
j£- > y. The growth of the instability is limited by the 
condition that u xx (x,h) > for all x at the free surface 
(the strain in the plane z = does not lead to any com- 
pression in the plane z — h), so that 9 — auk sin (uih) = 0. 
Thus, the amplitude of the vertical displacement, ka, is 
predicted to be proportional to 9 h. 

Discussion - Interestingly, the theoretical analysis, 
which only involves a decrease of the tensile stress associ- 
ated with the stretching of the material and an energetic 
cost associated with the induced shear, predicts that a 
stretched layer is always unstable. In agreement with 
the experimental observations, the instability does not 
exhibit any finite threshold, the amplitude of the modu- 
lation increases linearly with 9 for small strain (Fig. |4ji) 
and the wavelength A is proportional to the layer thick- 
ness h (Fig. [5]). In addition, the amplitude is predicted 
to be proportional to Oh, which explains why the pattern 
is more easily observed when the layer is thick and the 
stretching is large. 

For larger 9, the instability develops further and the 
deformation is no more sinusoidal. The stretching con- 
centrates in small regions in which we estimate that the 
relative displacement of two grains is of the order of 9X. 
When OX ~ 6, the typical size of a bridge, the material is 
locally torn off. Then, instead of an undulation of the free 
surface, one rather observes fractures. Thus, taking into 
account that A tx h and that 5 oc d, one predicts that 
thick layers or small grains, which are associated with 



very small critical strain ~ d/X, always lead to fractures, 
in agreement with the experimental observations. 

Finally, we note that the effect of the humidity content 
on the wavelength is accounted for by the dependence of 
A on the ratio £ = — (da/du xx )/G. For instance, in the 
limit of large u (small bridges), A ~ Ayj2a s /G0 m h. The 
increase of A with Rh would impose, in the framework of 
the crude model, that G increases slower than the ratio 
o's/Om = —(da/du xx ) (Note that <r s and m both increase 
with Rh)- However, the peculiar choice of a linear elastic 
response of the material to shear is certainly a matter 
of debate and we do no discuss this point further. Wc 
nevertheless point out that the most important features 
are not model dependent: the instability has in practice 
no threshold, the wavelength scales like the thickness and 
the amplitude of the undulation increases linearly with 
both the stretching and the layer thickness, . 

Conclusions - We reported for the first time the desta- 
bilization process of a material deposited onto an elastic 
membrane which exhibits a negative stress-strain sensi- 
tivity, da/du xx < 0, or equivalently "strain-softening", 
which generally leads to the fracturing of the material. 
However, as we demonstrated, the process differs ^signif- 
icantly from the usual fracturing processes 



ii, m 



as 

the instability develops gradually with increasing exter- 
nal strain. Interestingly, the instability mechanism re- 
ported here applies to any system having a negative sen- 
sitivity to stretching and, thus, is not expected to be 
specific to the granular matter. 
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